Pseudorapidity dependence of parton energy loss in relativistic heavy ion collisions 
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We analyze the recent data from the BRAHMS Collaboration on the pseudorapidity dependence 
of nuclear modification factors in Au+Au collisions at y^sjvjv = 200 GeV by using the full three di- 
mensional hydrodynamic simulations for the density effects on parton energy loss. We first compute 
the transverse spectra at t; = and 2.2, and next take a ratio Rr^ — RAA{r] = 2.2) / RAA{r] = 0), 
where Raa is a nuclear modification factor. It is shown that hydrodynamic components account 
for Rr^ ~ 1 at low px and that quenched pQCD components lead i?,, < 1 at high pT which are 
consistent with the data. Strong suppression at 77 
in the final state. 



2.2 is compatible with the parton energy loss 



PACS numbers: 24.85.+p,25.75.-q, 24.10.Nz 



Recent data from Relativistic Heavy Ion Collider 
(RHIC) reveal that hadron spectra at high px in cen- 
tral Au-t-Au colhsions are strongly suppressed relative to 
the scaled pp or large centrality spectra by the number 
of binary collisions [ij, |2, ylM l5| contrary to the enhance- 
ment in d-f-Au collisions [S H H 0- The RHIC data 
is consistent with the early predictions on jet quenching 
due to gluon bremsstrahlung induced by multiple scat- 
tering [lO| as a possible signal of dcconfincd nuclear mat- 
ter, the quark gluon plasma (QGP) (for a recent 
review, see Ref. 0|). Observed suppression of away- 
side peak in dihadron spectra in central Au-|-Au colli- 
sions 01 is also considered to be due to jet quenching, 
while correlation spectra in d-|-Au collisions are the same 
as in pp collisions [3 , where it is not expected to create 
hot and dense matter. Large elliptic flow observed in 
high-pT region is also considered to be a consequence of 
jet quenching 0,^3- Phenomenolpgical studies based 
on the parton energy loss llM ll^ are successful 

in describing various high px hadron spectra at RHIC: 
suppression of single particle spectr a V2^ fill li^ ^2^ . 
suppression of away-side correlation |23, |2J| , azimuthal 
anisotropy of high-pj hadron spectra in non-central 
collisions [2^ EM Eg Erf including centrality depen- 
dences [2^ . 

In addition to those data, the BRAHMS Collabora- 
tion recently reports the pseudorapidity dependence of 
the nuclear modification factors and shows that the yields 
of high pt charged hadrons are strongly suppressed even 
at 77 = 2.2 Furthermore, it is also shown that the 
ratio i?^^ = Rcp[ri = 2.2)/Rcp{'n = 0), where Rcp is 
a ratio of central to most peripheral yields normalized 
by the number of binary collisions, is almost unity at 
Pt < 2 GeV/c and R^^ < 1 at high px- These data 
are the first results of high px spectra in the forward 
rapidity region at RHIC and, thus, provide the novel op- 
portunity to study how the dense matter is distributed 
in the longitudinal directions. Therefore, further system- 
atic studies are necessary to confirm the presence of the 



jet quenching in the dense medium at RHIC. In this let- 
ter, we analyze p x sp ectra at p = 2.2 by employing the 
hydro-l-jet model [22, |2J, |23, |23| and test if the scenario 
of jet quenching in the QGP phase is still consistent with 
data. 

Hydrodynamics is found to be successful for the de- 
scription of the soft part of the matter produced in 
Au-|-Au collisions at RHIC especially in midrapidity re- 
gion {Y ~ 0) [23 ■ Motivated by these results, we also 
describe the space-time evolution of thermalized mat- 
ter even in off-midrapidity region {Y ^ 0) by solv- 
ing the equations for energy-momentum conservation 
in the full three-dimensional (3D) Bjorken coordinate 
(t, X, y, rjs) I^^jSI- Here r = y/t^ — z-^ is the proper time 
and ps = (1/2) ln[(t + z)/{t~ z)] is the space-time rapid- 
ity. Even at RHIC energies, one cannot observe "central 
plateau [s^l" in the rapidity distribution in Au-|-Au 
collisions. Note that a plateau-like structure in the pseu- 
dorapidity distribution observed at RHIC |3i US simply 
comes from the Jacobian between rapidity and pseudo- 
rapidity. In addition, elliptic flow as a function of pseu- 
dorapidity shows a peak at midrapidity [3^ . These data 
suggest that the full 3D hydrodynamic simulations are 
necessary for discussion on the global behavior in heavy 
ion collisions even at RHIC energies. Although our ini- 
tial condition for longitudinal flow rapidity (defined as 
Yf — (1/2) ln[(l -I- Vz)/{1 — Vz)] where Vz is the longitudi- 
nal flow velocity) is still assumed to be the exact scaling 
solution Yf{TQ, x,y,r]s) = Ps, we obtain the rapidity de- 
pendent observables by taking an p^ dependent initial 
energy density distribution which is factorized by a func- 
tion 



i?(ps) = exp 



(I r;. I -pflat/2)^ 



% I -?7flat/2) 



(1) 



where pflat and pcauss control the size of a flat (Bjorken- 
like) region near midrapidity and the width of Gaus- 
sian function in forward/backward rapidity respec- 
tively |37l |. The pseudorapidity distributions of charged 
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FIG. 1: (Color online) Invariant spectra of charged hadrons 
from UAl data IZlll in pp collisions at ^/s — 200 GeV is com- 
pared to PYTHIA predictions with Lund fragmentation (his- 
togram) and with independent fragmentation model (line), rj 
dependence is also shown for both models. Circles, squares, 
and triangles correspond to the predictions from PYTHIA 
with Lund fragmentation at = 0, 2.2 and 3.25, respectively 
which are compared to the results from PYTHIA with inde- 
pendent fragmentation model (lines). 



hadrons in central and semi-central collisions observed 
by BRAHMS |3| are already reproduced by choosing 
??flat/2 = 2 and yyoauss = 0.8 and also by taking account of 
local rapidity shift at each transverse coordinate 77so(a;, y) 
in Eq. |23. For details on how to parametrize ini- 
tial conditions in our hydrodynamic model, see Ref . |3ll | . 
In hydrodynamic calculations, a partial chemical equilib- 
rium (PCE) model with chemical freeze-out temperature 
T'^^ = 170 MeV is employed for the hadronic phase to de- 
scribe the early chemical freeze-out picture of hadronic 
matter [sJl. On the other hand, the QGP phase is as- 
sumed to be massless free partonic gas with the number 
of flavor Nf = 3. 

For the hard part of the model, we generate hard 
partons according to a pQCD parton model. We use 
PYTHIA 6.2 for the generation of momentum spec- 
trum of jets through 2^2 QCD hard processes. Initial 
and final state radiations are used to take into account 
the enhancement of higher-order contributions associated 
with multiple small-angle parton emission. The CTEQ5 
leading order parton distribution function |3fll | is used. 
Hadrons are obtained from an independent fragmenta- 
tion model option in PYTHIA. The ii'-factor K = 2.5, 
the scale Q = pT,jct/2 in the parton distribution function. 



and the primordial transverse momentum {k^)NN = 1-2 
GeV^/c^ are used to fit the neutral pion transverse spec- 
trum in pp colhsions at RHIC ^3- UAl data for 
charged hadrons are also well reproduced with this pa- 
rameters in the transverse momentum range pt > 2 
GeV/c as shown in Fig. ^ Independent fragmentation 
model is not applicable for the low pt transverse mo- 
mentum range. Hence we use the Lund string model for 
the pp reference when we compute nuclear modification 
factors at low pT (pt < 2 GeV/c). EKS98 parametriza- 
tion is employed to take into account nuclear shad- 
owing effect. The Cronin enhancement is modeled by the 
multiple initial state scatterings as in Ref. [i^ . 

Initial transverse positions of jets are determined ran- 
domly according to the number of binary collision dis- 
tribution. Initial longitudinal position of a parton is ap- 
proximated by the boost invariant distribution 13^ . Jets 
are freely propagated up to the initial time to(=0.6 fm/c) 
of hydrodynamic simulations by neglecting the possible 
interactions in the pre-thcrmalization stages. Jets are 
assumed to travel with straight line trajectory in the 
medium. 

Jets can suffer interaction with fluids and lose their 
energies. We employ the approximate first order formula 
(GLV formula) in opacity expansion from the reaction 
operator approach [l9j for the energy loss of partons in 
this work. The approximate first order formula in this 
approach can be written as 

^E = C drp (r, X (r)) (r - To) In (^^) . (2) 

Here we take L = 3 fm and /i = 0.5 GeV. C = 0.45 is an 
adjustable parameter and p(r, x) is a thermalized parton 
density in the local rest frame of fluid elements in the 
hydro-|-jet approach x (r) and Eq are the position 
and the initial energy of a jet, respectively. The initial 
energy Eq in Eq. ^ is Lorentz-boosted by the flow ve- 
locity and replaced by PqU^ where Pq and are the 
initial four momentum of a jet and a local fluid velocity 
respectively. The parameters related to the propagation 
of partons are obtained by fltting the nuclear modifica- 
tion factor for the neutral pion by PHENIX and are 
found to be consistent l24l with the back-to-back corre- 
lation data from STAR jia|. 

Let us start with the study of the space-time rapid- 
ity dependence of transverse dynamics from the full 3D 
hydrodynamics. In order to understand the dynamical 
effects on parton energy loss, we plot in Fig. |21 parton 
densities and 7iong-weighted transverse flow velocities av- 
eraged over the QGP and mixed phases as a function of 
T at = 0, 2, and 3. Here a gamma factor of longi- 
tudinal flow is 7iong = cosh(yf). Note that the quan- 
tity (wT7iong) is independent of the space-time rapidity 
77s when the initial longitudinal flow and the initial en- 
ergy density obey the scaling solution If (tq, x, y, 77s) = 77s 
and e(ro, x, y, r]s) = e(To, x, y) respectively. We note that 
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FIG. 2: (Color online) Time evolutions of average parton 
density (upper panel) and average 7iong-weighted transverse 
velocity (lower panel) at space-time rapidity r^s ~ 0, 2, and 
3 from hydrodynamic simulations in 0-10% central Au+Au 
collisions at y^ijvjv = 200 GeV. Impact parameter is 6 = 3.7 
fm. Average is taken only in the QGP and mixed phases. 



longitudinal flow from full 3D hydrodynamics remains to 
be close to the Bjorken flow; Yf « 77s at t^s < 3 within 
our initial conditions used in this work. As shown in the 
flgure, the difference of the time evolutions of the ther- 
malized parton density and the transverse flow velocity 
between 77s = and ?7s = 2 are practically tiny, since the 
initial scaling region reaches to ijs — ??flat/2 — 2 [isj . 
Thus, the dynamical effect on the jet energy loss at 
r] = and 2.2 is expected to be the same, whereas the 
amount of jet energy loss at 77 = 3 must be small be- 
cause of the smaller parton density. As will be shown in 
Fig. 21 we should emphasize here that this does not mean 

We next show in Fig.|21the transverse momentum dis- 
tributions for charged hadrons from the hydro-|-jet model 
in central Au-|-Au collisions at RHIC. Thermal freeze- 
out temperature T^^ = 100 MeV is used in the calcula- 
tion. This choice is consistent with the PHENIX data at 
,/snn = 130 GeV [sJl. Each spectrum is the sum of the 
soft component and the hard component. Before sum- 
mation, the hard component is multiplied by a "switch" 
function [i^ {1 -|- tanh[2(pT — 2)]}/2 (where pr is in the 
unit of GeV /c) in order to cut the unreliable components 
from the independent fragmentation scheme and also to 
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FIG. 3: (Color online) Transverse momentum distributions 
of charged hadrons in Au -I- Au collisions at -^/ijvjv = 200 GeV 
is compared to data from BRAHMS 0. Solid lines represent 
the hydro-|-jet result averaged over I77I < 0.1, 2.1 < I77I < 
2.3, and 3.0 < |r;| < 3.5. The impact parameter for 0-10% 
centrality is taken to be fe = 3.7 fm. In the bottom panel, the 
ratio of the spectrum at 77 = 2.2 to that of 7; = is plotted. 
Enhancement at very low pt comes from the effect of Jacobian 
for the transformation of rapidity to pseudorapidity. 



obtain the smooth spectra. The hydrodynamic compo- 
nents are dominated in the range of pt < 2 GeV/c. The 
slope of hadrons in low px region at ry = is nearly the 
same as the one at 77 = 2.2 as clearly seen in the bot- 
tom panel of Fig. |3| in which the ratio of the spectrum 
at 77 = 2.2 to that at 77 = is plotted. This results from 
the similarity of transverse dynamics between 77 = and 
2 as shown in Fig. |21 On the other hand, the slope from 
pQCD components in high region becomes steeper as 
77 increases reflecting the original pp spectra (see Fig. ^ . 

We now turn to the study of the nuclear modification 
factors Raa for charged hadrons defined by 



dN- 



Raa 



AT dN'P + 'P ' 

J V coll rf2p^^^ 



(3) 



where A^coii is a number of binary collisions. Figure 0] 
shows the nuclear modification factors Raa for charged 
hadrons at 77 = 0, 2.2, and 3.25 in Au+Au collisions at 
RHIC for an impact parameter b = 3.7 fm. The nu- 
clear modification factor Raa^s in low px region (px/S 2 
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FIG. 4; (Color online) Nuclear modification factors are 
compared to the BRAHMS data in Au+Au collisions at 
= 200 GeV. 
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FIG. 5: (Color online) Elliptic flow parameter V2 for charged 
hadrons at 77 = 0, 2.2, and 3.25. The impact parameter is 
taken as 6 = 7.2 fm. We only take into account contributions 
from hard components. 



GcV/c). where the hydrodynamic component dominates, 
at 7^ = and 2.2 are almost identical. This is due to the 
comparable time evolution of the parton density at = 
and 2.2 in hydrodynamics as shown in Fig. [21 PYTHIA 
prediction reveals that the spectrum at low pT is very 
similar within 77 < 2 up to a factor of 30% in pp colli- 
sions as shown in Fig.^ Raa{v = 0) > Raa{ii = 2.2) at 
high pt is a consequence of the steeper slope at 77 = 2.2 
compared to the slope at 77 = in pQCD. When the pr 
slope is steep, the nuclear modification factor becomes 
sensitive to nuclear effects: a small shift of a spectrum is 
likely to produce a large effect on the ratio of the shifted 
spectrum to the original one. It should be noted that, 
due to the above reason, the nuclear modification fac- 
tor from the Cronin effect at SPS energies j4g is much 
larger than the one at RHIC energies 013, SIJ- The nu- 
clear modification factor at 7y = 3.25 in the range < 5 
GeV/c is larger than at midrapidity, because thermal- 
ized parton density at 7; = 3.25 is about 40% smaller 
than at midrapidity. However, RAAiv = 3.25) eventually 
becomes smaller than the one for 77 = or 2.2 in high pT 
region. This is due to the much steeper slope at high px- 
We also plot the ratio i?^ defined by 



RAAiv = 2.2) 
RAAiv = 0-Oy 



(4) 



which can be compared with the data from BRAHMS 
defined by R^^ = Rcp{v = 2.2)/Rcp{ti = 0.0). Hy- 
drodynamic predictions at very peripheral collisions are 
likely to be unreliable, since parton density is too small 
to justify the assumption of the local thermal system in 
Au-I-Au collisions at RHIC. In Fig.^ i?^ is found to be 
almost unity in low region and gradually decreases 
with px as consistent with data as expected in the above 
consideration. 



One may worry about the nuclear shadowing effect on 
the since the accessible range of the parton momen- 
tum fraction x in the parton distribution functions at 
77 = 2.2 is smaller than at midrapidity. However, nu- 
clear shadowing effect on Raa{v = 2.2) < RAAiv = 0) 
is found to be ~ 15 % around px = 2-3 GeV/c within 
EKS98 parameterization. Other models for nuclear shad- 
owing should be investigated elsewhere. 

In order to disentangle the density effect from the 
"slope" effect among nuclear modification factors, we 
predict the elliptic flow parameter V2 = (cos(2(/))) for 
high px charged hadrons in forward rapidity region as 
shown in Fig. |31 In this calculation, we only take ac- 
count of the contributions from hard components. We 
choose the impact parameter as = 7.2 fm correspond- 
ing to 20-30% ccntrality. V2 in high px region is gener- 
ated by jet quenching [25|. The difference of path length 
causes the difference of the amount of parton energy loss 
in azimuthal directions in non-central collisions. Con- 
sequently, one can expect from Eq. ||2J) that the higher 
parton density results in the larger positive V2 until V2 
reaches the limiting value 01 ■ we expected, our pre- 
diction 112(7; = 3.25) < V2{r] = 0) is a clear evidence of the 
density effect on the parton energy loss. We comment on 
our prediction on the V2 at high px > 5 GeV/c at midra- 
pidity. Our prediction is similar to the results of the 
recent work |48| in which they found that the azimuthal 
anisotropy of high px particles underestimate for a real- 
istic nuclear density profile although hard-sphere nuclear 
profile looks consistent with the data [23l |. 

In summary, we have studied the pseudorapidity de- 
pendence of the nuclear modification factors for charged 
hadrons within the hydro-|-jet model. As well as the yield 
of charged hadrons, the radial fiow at t/s = 2 is found to 
be very similar to that at 77s = 0. This results in i?,; ~ 1 



5 



at pt < 2 GeV/c. i?^ < 1 in high region can be un- 
derstood by the steeper slope of pr spectrum at -q — 2.2 
than at ?7 = from the pQCD components. This suggests 
that the longitudinal region of dense partonic matter pro- 
duced in Au-|-Au collisions reaches to ?7s ~ 2 and that 
strong hadron suppression at off-midrapidity is consistent 
with the final state parton energy loss in the medium. 
This reminds us the previous analysis of pseudor apid ity 
dependence of elliptic flow at the RHIC energy [13,01 in 
which elliptic flow can be reproduced by hydrodynamics 
only in the region | 77 | ^ 2. We also predicted the elliptic 
flow parameter V2 for high pT charged hadrons in both 
midrapidity and forward rapidity regions. We found the 
strong rapidity dependence of V2{pt)- 

It would be very interesting to see how robust the 
present results arc in more realistic analysis. For exam- 
ple, the lattice studies predict the non-trivial behavior 
of screening mass in the vicinity of phase transition re- 
gion This affects the space-time evolution of the 
parton density in the energy loss formula in Eq. Q . The 
parton density is multiplied by the correction factor for 
the effective degree of freedom ^^//ig |13 where and 
/i are, respectively, perturbative and non-perturbative 
screening mass. The reduction of energy loss due to the 
non-perturbative behavior of the screening mass is ex- 
pected to be large in forward rapidity region where the 
initial parton density is not so large. This modification 
causes the change of equation of state (EOS) as well as 
the amount of the parton energy loss. Hence we may 
need to use different initial conditions and simulate the 
hydrodynamic model with a more realistic EOS. This is 
beyond the scope of this paper although it should be re- 
visited elsewhere. 
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FG03-93ER40792. 
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